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Abstract 

In this paper, we study dynamics of the charged plane symmetric 
gravitational collapse. For this purpose, we discuss non-adiabatic flow 
of a viscous fluid and deduce the results for adiabatic case. The Ein- 
stein and Maxwell field equations are formulated for general plane 
symmetric spacetime in the interior. Junction conditions between 
the interior and exterior regions are derived. For the non-adiabatic 
case, the exterior is taken as plane symmetric charged Vaidya space- 
time while for the adiabatic case, it is described by plane Reissner- 
Nordstrom spacetime. Using Misner and Sharp formalism, we obtain 
dynamical equations to investigate the effects of different forces over 
the rate of collapse. In non-adiabatic dynamical equation is 

joined with transport equation of heat flux. Finally, a relation between 
the Weyl tensor and energy density is found. 

Keywords: Gravitational collapse; Junction conditions; Dynamical equa- 
tions; Transport equations. 



1 Introduction 

Gravitational collapse is defined as the implosion of a celestial body under the 
influence of its own gravity. It is one of the basic processes driving evolution 
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within galaxies, assembling giant molecular clouds and producing stars. The 
study of gravitational collapse is motivated by the fact that it represents one 
of the few observable phenomena in the universe. The gravitational collapse 
of a star proceeds to form compact objects as white dwarf, neutron star or 
black hole primarily depending upon the mass of the star. The end state 
of continual gravitational collapse of a massive star is an important issue in 
gravitation theory. According to Penrose, gravitational collapse of a massive 
star gives rise to a spacetime singularity [T]. 

Oppenheimer and Snyder [2] innovated the first mathematical model for 
the description of gravitational collapse of stars. They dealt with the dust 
case and concluded that gravitational collapse of massive stars might end to 
form a black hole. This work of Oppenheimer and Snyder attracted many 
people to explore it further. Markovic and Shapiro [3] extrapolated pio- 
neer's work involving positive cosmological constant and explained its effects 
on the rate of collapse. Misner and Sharp [I] did the same job by taking 
ideal fluid and also discussed some other aspects like thermodynamics and 
hydrodynamics of gravitational collapse. 

The standard approach dealing with gravitational collapse problems re- 
quires junction conditions to join the two spacetimes. Three types of junction 
conditions are available in literature proposed by Darmois [5], Lichnerowicz 
[6] and O'Brien and Synge [?]. Bonnor and Vickers [8] studied these junc- 
tion conditions and showed that the conditions suggested by Darmois and 
Lichnerowicz are equivalent. Also, the O'Brien and Synge conditions are un- 
satisfactory. It was concluded that Darmois junction conditions are the most 
convenient and appropriate. 

Most of the work in this scenario has been done using spherical symmetric 
models [9]-[TT]. Herrera et al. [12] used complicated fluid for this purpose and 
discussed the physical meaning of expansion free fluid evolution. In another 
paper [13] . Herrera et al. explored the dynamics of dissipative gravitational 
collapse considering spherical symmetry. They also discussed applications of 
the results to astrophysical scenario. 

The phenomenon of gravitational collapse has also been explored in the 
context of cylindrical, planar and quasi-spherical symmetries [2]-[!l9|. Kurita 
and Nakao [20] studied the null dust collapse for the cylindrically symmetric 
spacetimes. They concluded that singularity must form at symmetry axis and 
also discussed the behavior of geodesies arriving at the singularity. Wang et 
al. [21] worked on self-similar plane symmetric solutions. It was found that 
during collapse trapped surfaces are not formed and singularities are space- 
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like. Nath et al. [22] elaborated junction conditions for the quasi-spherical 
Szekeres spacetime in the interior and the Riessner Nordstrom Vaidya space- 
time in the exterior region. The time difference between the formation of 
apparent horizon and central singularity was also discussed. 

Using the concept of gravitational lensing (GL), Virbhadra et al. [23] in- 
troduced a new tool for examining naked singularities. Gravitational lensing 
is the process of bending of light around a massive object such as a black hole. 
Virbhadra and Ellis [23] discussed GL by the Schwarzschild black hole. It was 
found that the relativistic images guarantee the Schwarzschild geometry close 
to event horizon. The same authors [25] also analyzed GL by a naked singu- 
larity. Claudel et al. [2S] proved that the necessary and sufficient condition 
for the black hole to be surrounded by a photon sphere is that a reasonable 
energy condition holds. Virbhadra and Keeton [27] showed that weak cosmic 
censorship hypothesis (CCH) can be examined observationally without any 
uncertainty. Virbhadra [21] found that Seifert's conjecture is supported by 
the naked singularities forming during Vaidya null dust collapse. The same 
author developed an improved form of CCH using GL phenomenon [29J 

In a recent paper, Sharif and Abbas [5U] found the effects of electromag- 
netic field on the gravitational collapse for perfect fluid in the presence of 
cosmological constant. It was concluded that charge increases the rate of 
collapse by decreasing the limit of cosmological constant. The same authors 
extended this work for 5D collapse [3T] and found that the range of appar- 
ent horizon is greater than 4D case. Di Prisco et al. [32] figured out the 
consequences of charge and dissipation for spherical symmetric gravitational 
collapse of a real fluid. It was assumed that heat flow, free streaming ra- 
diation and shearing viscosity are the causes of dissipation. The dynamical 
and transport equations are coupled to observe the effects of dissipation over 
collapsing process. 

This paper extends the work of Di Prisco et al. [32] to plane symmetry. 
The models exhibiting plane symmetry may swear out as test-bed for nu- 
merical relativity, quantum gravity and contribute for examining CCH and 
hoop conjecture among other important issues. The paper is organized in the 
following pattern. In the next section, the dynamics for non-adiabatic flow 
is discussed. The Einstein-Maxwell field equations and junction conditions 
are found. The dynamical and transport equations are obtained and then 
coupled. A relation between energy homogeneity and the Weyl tensor is also 
given. Section 3 describes dynamical equations for adiabatic case. The last 
section summarizes the results. 
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2 Dynamics of Viscous Non-adiabatic Case 



Here the dynamics for non-adiabatic flow is discussed. We formulate dynam- 
ical and transport equations and then finally couple these equations. 

2.1 Interior Spacetime and Some Physical Quantities 

We consider a plane symmetric distribution of collapsing fluid bounded by 
a hypersurface S. The line element for the interior region has the following 
form 

ds 2 _ = -A 2 (t, z)dt 2 + B 2 (t : z)(dx 2 + dy 2 ) + C 2 (t, z)dz 2 : , (2.1) 

where we have assumed co- moving coordinates inside S. The interior coor- 
dinates are taken as x~° — t> X~ l — x -, X~ 2 — X~ 3 — z - It is assumed 
that fluid is locally anisotropic and suffering dissipation in the form of shear- 
ing viscosity, heat flow and free streaming radiation. The energy-momentum 
tensor has the following form 

T ab = {fJL + P±)V a V b + P ± g ob + (P g -P ± )xaXb + q a Vb + V a q b + ee a e b 
- 2v<r ab , (2.2) 

where fx is the energy density, P z is the pressure in ^-direction, P± is the 
pressure perpendicular to z-direction (i.e., x or y direction), V a four velocity 
of fluid, r) is the coefficient of shear viscosity, \ a is a unit vector in z-direction, 
q a is the heat flux, e is the radiation density and i a is a null four vector in 
z-direction. Furthermore, these quantities satisfy the relations 

V a V a = -1, X a Xa = l, X a V a = 0, 

V a q a = 0, £ a V a = -l, £ a i a = 0. (2.3) 

Since the metric defined in Eq. fl2.ip is co-moving, we can take 

V a = A~ l 5l X a = C- x 5l q a = q C-Hl 
t = A- l 5 a + C~ l 5l. (2.4) 

In the standard irreversible thermodynamics by Eckart, we have the fol- 
lowing relation [S3] 

7t ab = -2 Wab , n = -ce, (2.5) 

where rj and ( stand for coefficients of shear and bulk viscosity, a ab is the 
shear tensor, G is the expansion, IT is the bulk viscosity and ir ab is the shear 
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viscosity tensor. The algebraic nature of Eckart constitutive equations causes 
several problems but we are concerned with the causal approach of dissipative 
variables. Thus we would not assume (12. 5p rather we shall resort to transport 
equations of Miiller- Israel-Stewart theory. 

The non-zero component of acceleration and expansion scalar are given 

by 

A' „ 1,25 C\ 

a: > = T b= a { 1T + c ) ' < 2 ' 6 > 

where dot and prime denote differentiation with respect to t and z respec- 
tively. The non-vanishing components of shear tensor are 

*n = ~B 2 F = a 22 , a 33 = ^C 2 F; F=j(^- + ^). (2-7) 
The magnitude of the shear tensor, i.e., the shear scalar a is defined as 

* 2 = \o*°* = \f 2 (2.8) 

which implies that F 2 = 9a 2 . 

2.2 The Einstein and Maxwell Field Equations 

The energy-momentum tensor of an electromagnetic field is defined as 

E ab = ^{F c a F bc - \F cd F cd g ab ), (2.9) 
where F ab is the electromagnetic field tensor given by 

F ab = <Pb,a - <Pa,b, (2.10) 

(f> a is four potential. The Maxwell field equations are 

F ab ;b = fi r, F [ab , c] =0, (2.11) 

where fiQ = An is the magnetic permeability and J a is the four current. In a 
co-moving frame (p a and J a are given by 
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where £ and </> are the charge density and electric scalar potential respectively 
and both are functions of t and z. The charge conservation J a . a = gives 
the charge for interior region as 

s{z) = [ iB 2 Cdz. (2.13) 



^AC\ (2.14) 



For a = 0, 3, the first of Eq. fl2.lip becomes 
^ [ A C B w 



+ = °' (2 - 15) 

while for a = 1,2, it is trivially satisfied. Also, the second of Eq. fl2.lip is 
identically satisfied. Integration of Eq.f l2.14p with respect to z, assuming 
(f>'(t,Q) = 0, gives 

fips(z)AC 

4> = — 2 • (2.16) 

The Taub's mass function [34J for plane symmetric spacetime can be gener- 
alized to include the electromagnetic contributions as 

B , B 2 B'^ s 2 , 

™(M) = - or) + w (2.17) 

For the interior spacetime, the Einstein field equations, G a b = 87r(T a {, + 
E a b), yield the following set of equations 

^ + «K + ^! = |(^ + | ) + ( ^[^: 

+(£ -f)f], P.18) 

i?' j4'S CB' 
-8^C(, + £ ) = -- + — + — , (2.19) 

BMP, + ? ,F)^ + (*|-)> = -(-)*[- + 5-^ 

C, BC, B 2 .A" B" A'C B\ B'C\ .„ „„. 
+ C )+ BC l + C? l T + lf" ( C-B ) - BC 1 ' (2 - 20) 

Mft + « - - ^ = -(^) 2 [f + 

2AB. 



AB B AB 
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In view of Eqs.f l2.6j) and ( 12. 7p . Eq.f l2.19j) can be written as 

47r(g + e)C=i(e-*y-F^. (2.22) 

2.3 Junction Conditions 

Here we formulate junction conditions for the general plane symmetric space- 
time in the interior and charged Vaidya plane symmetric spacetime in the 
exterior. The line element for the exterior region is [35] 

dsl = (^MM. _ ( -^)dv 2 - IdZdv + Z 2 {dX 2 + dY 2 ), (2.23) 

z z 

where x + ° = v i X +1 = %■■> X +2 = ^ X +3 = The metric for hypersurface 
is defined as 

(ds% = ~{drf + f{rf{dx 2 + dy 2 ), (2.24) 

where w % = (r,x,y)(i = 0,1,2) are the intrinsic coordinates of S. The 
equations of hypersurface in terms of interior and exterior coordinates are 

kr(t,z) = z-z s = 0, (2.25) 
k + (v,Z) = Z-Zx(v) = 0, (2.26) 

where is constant. Using Eqs. 02.25]) and ( I2.26p . we get the interior and 
exterior metrics over the hypersurface as 



(ds 2 _)j; = -A 2 {t,zx)dt 2 + B 2 {t,zx){dx 2 + dy 2 ), (2.27) 

(dslh = _ [( =^) + !j)) + ^ 

+ Zl(dX 2 + dY 2 ). (2.28) 

Now we use the junction conditions proposed by Darmois [5], the first 
condition is 

(ds% = (ds 2 _h = (ds 2 + h (2.29) 
which yields the following equations 

dt 1 
d^ ~ A' 



(2.30) 



Z s = B, (2.31) 
du .,-2Mf» e 2 (u). 2dZ : 



[(^P^ + ^Pr) + ==7^ • (2-32) 



dr Zy, Z\ dv 



The second junction condition is the continuity of extrinsic curvature (the 
second fundamental form) 



K ab = K~ b = K+. (2.33) 

The unit normal in terms of interior and exterior coordinates are given re- 
spectively as 

n~ = C(0, 0,0,1), n + a = (-i S) 0, 0, P), (2.34) 

here hat denotes differentiation with respect to r. The surviving components 
of the extrinsic curvature for the interior spacetime are 

K 00 = -[^h, K' x = [^] s = K-. (2.35) 

The non-null components of the extrinsic curvature for the exterior spacetime 
are given by 

,d 2 v ,dv . i ,M e 1 . ,dv ^ 

A » = teW 1 ~ ~ < 2 - 36 > 

K+ = [Z^- 2Af)^] E = A'+. (2.37) 

Thus the second junction condition yields the following equations 

< A ' h = [£(*)->_(".£)(£)],,, (2.38) 



L AC> y dT 2y dr J V Z 2 Z^'^dr' 

r BB' . dZ ,e 2 .dv, 

I— 1= = l z * + <z- 2M >*]- < 2 ' 39) 

After some algebra, it follows that 

M(v)=m(t,z) ^ s = e, (2.40) 

9 = ^-^-^(^0-1)- (2-41) 

These equations give necessary and sufficient conditions for the matching of 
interior and exterior spacetimes. Equation ( I2.4ip describes a relationship 
between heat flux, effective pressure in ^-direction and charge over the hy- 
persurface. It shows that if the fluid has no charge then effective pressure 
and heat flux are equal over the hypersurface. 
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2.4 Dynamical Equations 

Here we develop equations that govern the dynamics of non-adiabatic plane 
symmetric collapsing process by using Misner-Sharp formalism |4]. The 
proper time derivative and proper derivative in ^-direction are defined re- 
spectively [IS] as 

D * = j>i' (2 - 43) 

where Z = B. The velocity of the collapsing fluid is the proper time deriva- 
tive of Z defined as 

U = D f (Z) = | (2.44) 
which is always negative. Using this expression, Eq. ()2.17p implies that 

B' r ~ 2m s 2 ,i , 

E -c = [u "ir + i^ < 2 - 45 > 

When we make use of Eq. (I2.43[) in (I2.22p , we have 

4nC(q + e) = E[±D 2 (G - F) - I]. (2.46) 

The rate of change of mass (given in Eq.f l2.17f ) with respect to proper time 
is given by 

Df m = -An[(P z + e - ^ V F)U + E(q + e)}Z 2 + g(^ _ i). ( 2 .47) 

This equation shows how mass is varying within the plane hypersurface 
under the influence of matter variables. The first term represents effective 
pressure in z-direction and radiation density. When collapse takes place, 
this term is positive implying that energy increases by this factor. The 
second term in square brackets shows that energy is going out from the 
plane hypersurface while the last term is the charge contribution. During 
collapse, energy decreases due to these terms. Similarly, we calculate 

D^m = 47r[// + e + (g + e )|]Z 2 + Id 2 s + i(^ - 1). (2.48) 
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This equation describes how different quantities influence the mass between 
neighboring hypersurfaces in the fluid distribution. The term (fi+e) indicates 
the effects of energy density and radiation density. Similarly, the second term 
shows the amount of heat and radiation which is getting out. The remaining 
two terms represent contribution of electric charge. Integration of Eq. fl2.48p 
yields 

r 2 u ~ s 2 1 r 2 s 2 

771 = J Q ^ + t + (l + ^ Z2dZ+ ^ + 2j { J^ )dZ 

+ { -^^j\t)d~Z. (2.49) 

The dynamical equations can be obtained from the contracted Bianchi 
identities (T ab + E ab ).b = 0. Consider the following two equations 

(T ab + E ab ), b V a = (T ob + Ef)V = 0, (2.50) 
(T ab + E ab ), bXa = (if + Ef b ) X3 = (2.51) 



which yield 



^±± + ( fl + 2 e + P,-^,Fl-^ + 2( lt + e + Pj. + ^ l F) 
* ■ W + ^( 5 + £ ) + ^( 9 + £ ) = 0, (2.52) 



AB AC AC y ^ ' BC 

4 A' . B' uhs' 

The acceleration of the collapsing fluid is defined as 

Using Eqs.(j22U), (TX42]) and (l2~T7j) . we have 

4 ~ m s 2 EA' 
D f U = -4tt(P z + e - - V F)Z - ^ + ^(/i 2 + 1) + — (2.55) 
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which gives the value of 4^ 



4 = ^[D f U + An(P z + e - -r)F)Z\ + ^S- - + 1). (2.56) 

A £ L T v 3' ; J £Z 2 2EZ 3 



Substituting this value in Eq. (12.53 p . it follows that 




This equation yields the effect of different forces on the collapsing process. 
It can be interpreted in the form of Newton's second law of motion i.e., force 
= mass density x acceleration. The term within round brackets on LHS 
represents the inertial or passive gravitational mass density. This term shows 
that effective pressure, energy density and density of null fluid have effects 
on mass density while heat flux and charge have no contribution here. By 
equivalence principle, the round brackets factor on RHS is taken as active 
gravitational mass density. The quantities within square brackets in the first 
term show the influence of effective pressure, radiation density and electric 
charge on active gravitational mass. Using Eq. fl2.49p in (I2.57p . it follows that 
charge increases the active gravitational mass if 



If this inequality holds then charge regeneration phenomenon analogous to 
pressure regeneration occurs [36]. The pressure regeneration means that the 
pressure which is trying to keep hydrostatic equilibrium by balancing gravi- 
tational attraction, at the same time contributes to the active gravitational 
mass. This implies that it promotes gravitational collapse. Otherwise, if 
the above inequality is not satisfied, charge will decrease active gravitational 
mass and consequently the Coulomb repulsion may prevent the gravitational 
collapse. 

The first term in the second square brackets is the gradient of effective 
pressure in ^-direction and radiation density. Since this gradient is negative, 




(2.58) 
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it increases the rate of collapse. The second term is due to local anisotropy 
of pressure, radiation density and contribution of viscosity. If this term is 
positive then it contributes to increase collapse and vice versa. The last 
term depicts Coulomb repulsion that opposes gravitation implying that it 
decelerates the collapsing process. 

Finally, the last square brackets is entirely due to dissipation. To see the 
role of Dfq, this equation is coupled with causal transport equation. The 
consequences of Dfe have been discussed by Misner [37]. The outward flux 
of radiation accelerates collapsing process by increasing gravitational force. 
The third term is positive as U < 0, so it slows down rate of collapse. The last 
term shows the combine effect of viscosity and dissipation. From Eq. fl2.57p . 
the condition for hydrostatic equilibrium can be obtained by replacing U = 
q = e = r] = 0as 

D- 7 P Z = - Z h =~ + 1 + ~ - MP Z - Pa.). 2.59 

2.5 Transport Equation 

The transport equation for heat flux derived from the Miiller-Israel-Stewart 
theory of dissipative fluids [33J is given by 



r h ab V c q b , c + q a = -Kh ab (T b + Ta b ) - ± K T 2 (^) ]b q a , (2.60) 

where h ab is the projection tensor, k denotes thermal conductivity, T is tem- 
perature and t stands for relaxation time. This equation has only one inde- 
pendent component 

nT 2 q to 3U F 1 kE kT 

D T q = -^^ D f(-^)- ( l\^ + ^ + -\-— D Z T -—^ 

2r ki z 2Z 2 r r r h 

kT 4 ~ s 2 1 

x D f U-—[m + A>ir(P z + e-- V F)Z 3 -^]j- 2 . (2.61) 

We now couple this equation with dynamical Eq. (12. 57|) to see the effects of 
heat flux or dissipation on collapsing process. Replacing Eq. fl2.6ip in (I2.57p . 
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we obtain 

kE 2 



(H + P z + 2e - -r]F){l - a)D f U = (1 - a)F grav + F hyd + 



D T + (JjL) _ D e ] _ + -F - — ) - 2Fe 

z 1 2r tv kT 2; t j yv 2Z 2 r ; 

(^ + n (2-62) 

where F grav , Fh y d and a are given by the following equations 

F g rav = -(Li + P z + 2e-^ V F)[m + 4ir(P z + e-^- V F)Z 3 -^L 

x (Aig + 1)]^, (2.63) 

F M = -F 2 [D^(F, + 2 e -^F) + |(P 2 -F ± + e-2r ? F)-^|J], 

(2.64) 

kT 4 
a = —^ + P z + 2e--r 1 F)-\ (2.65) 
r 3 

The consequence of coupling transport and dynamical equations is that the 
inertial and active gravitational mass densities are affected by a factor a 
given by Eq. (12.65 p . The gravitational force term defined in Eq. (12.631) is also 
affected by a but the hydrodynamical forces Eq. ()2.64p are not influenced by 
this term. 



2.6 Relation Between the Weyl Tensor and Matter 
Variables 

Here we find some relationship between the Weyl tensor and matter variables. 
The Weyl scalar C in terms of Kretchman scalar TZ, the Ricci tensor R ab and 
the Ricci scalar R is given by 

C 2 = K - 2R ab R ab + -Pi 2 . (2.66) 

3 

The Kretchman scalar TZ 2 = R abcd R a b cl i becomes 

2 1 1 

+ ^^ 2323 ) 2 -I^( i?0113 ) 2 ]- ( 2 - 67 ) 
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The non-zero components of the Riemann tensor can be written in terms of 
the Einstein tensor and mass function as 



1 Is 2 
R 0l0l = (AB) 2 i—G 33 + —(m-—)} = R 0202 , 

2 1 1 1 2 s 2 

#0303 = (AC) [7^ G 00 - ^2~ G 33 + J^ G 22 ~ ^ ( m ~ ^g)]> 

R 12l2 = 2B(m-^), 

#1313 = (#G) 2 [^G 00 - ±(m - ^)] = # 2323 , 

#0113 = — ^ — ^ 03 ' 
Substituting these values in Eq. fl2.67p . after some algebra, we obtain 



48 s 2 2 16 s 2 Goo G 33 G 



S 6V 2S y B 3 " 2B A 2 C 2 B 21 

4 n 2 , qr/^00x 2 , /C33x2i , ^^2 ^00^33 
^2(^2^03 + 6 li A2 ) c2 ) \ + B A U 22 ^ A2( j 2 

Goo G33 G 2 2 (o RQ\ 

+ 4 (^-^)^- ( 2 - 68 ) 

Now we calculate the remaining part of the Weyl scalars which need Ricci 
tensor and Ricci scalar in terms of the Einstein tensor. These are 

d _ A 2\ G oo , G 33 G 22 

~ 2A 2 2G 2 "ff 2 -"03-^03, 

r>2 /~i rt s~i r~i rt 

D Lr (^33n _ D r> _ ^r^OO , ^33 ^22-, 

■n-ii — Tri~r^ — ^22, -K33 — t> 



# 



2 A 2 C 2 ' 00 2A 2 2G 2 5 2 

Goo G33 2G22 
"A 2 " _ ~C 2 ~ B 2 ' 



^-y2 O/^ 2 9^ 2 

^>a6^> _ ^00 1 ^33 1 ^22 ~ lT n:; 



A 4 G 4 5 4 A 2 G 2 ' 

Thus the remaining part of the Weyl scalar becomes 

Iff 2 _ iR^R 5 Gq 5 G33 8 G 2 2 4Gp 3 2 G00G33 

3 3 A 4 3 G 4 3B 4 A 2 C 2 3 A 2 G 2 

4 G22G33 4 G00G22 /„ „„x 

+ s"^ 2- "^^ 2 ^ 2- - 1 j 
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Using Eqs.( l2.68"j) and (12.691) . the Weyl scalar takes the form 

n 2 48 s 2 2 16 s 2 G o G 33 G22, 



fi6 v 2 P/ S 3V 2B n A 2 C 2 B 21 

4 n 2 , or/^00\2 1 /^33\2] , jLr 12 o Gq G;; 
^2^2^03 + J K ^ 2 ^ + ^ (J2 I \ + ^4^22 Z ^ 2( ^ 

. /Goo _ G33. G22 _ 5 (jqq ^ 5 G§ 3 ^ 8 G 2 2 1 4G 2 

V /1 2 /^2 ' D2 o AA o nA o IDA 



A 2 C 2 B 2 3 A 4 3 G 4 3B 4 A 2 G 2 

2 G00G33 4 G22C33 ^ 4 G00G22 

3 A 2 G 2 3 C 2 B 2 ~ 3 A 2 S 2 ' 
After some algebra, it leads to the following equation 



(2.70) 



CB 3 s 2 £> 3 Goo . G 2 2 G33 3s 2 /r, ni\ 



(48)5 v S y 6 A 2 P 2 G 2 P 4 
Using the field equations, we have 
Goo . G 2 2 G33 3/1qS 2 

^42- + ^2---^2---^4- = 87r (/ i + jP4 --^ + 2r ' F )- ( 2 - 72 ) 

In view of the above equation and using Z = B, Eq. (12.711) becomes 

CZ 3 An ~ s 2 

[m-—(fi + P ± -P z + 2r]F)Z 3 --^^ 2 + l)]. (2.73) 



(48)5 3 ^ - 2Z 

The derivatives of ( cz \ ) with respect to T and Z are given by 

(48)2 

CZ 3 1 2 

^(7^7t) = -^kZ 3 D f ( f i + P ± -P z + 2 V F) + ( f i + P ± + e+-r ) F) 
(48)5 3 3 

x Z 2 U + ( q + e)EZ 2 } + S ^- 2 ( f 4 + l), (2.74) 

CZ 3 Z 2 U 1 

= 4vr[(g + e)^--Z 2 ^( / , + P ± -P 2 + 2r / F) + ( e -P ± 
(48)5 E 3 

+ P 2 - 2?? F)Z 2 ] - S -^f{^ + 1) + -t + 1). (2.75) 

Z 2Z A 

These equations give relationship between the Weyl scalar and the fluid prop- 
erties like density, viscosity and pressure (anisotropy). For perfect and non- 
charged fluid, Eq. (1 2.75ft reduces to the following form 

C 7 3 A.7T 

Dz(j^rr) = ^Z 3 D^. (2.76) 
(48)5 3 
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Using the regularity condition, it is concluded that -D^/i = if and only 
if C = 0. This means that if energy density is homogeneous, the metric is 
conformally flat and vice versa. 

We would like to mention here that a particularly simple relation between 
the Weyl tensor and density inhomogeneity such as (2.76), for perfect non- 
charged fluids, is at the origin of Penrose's proposal to provide a gravitational 
arrow of time in terms of the Weyl tensor j3B] • The rationale behind this idea 
is that tidal forces tend to make the gravitating fluid more inhomogeneous 
as the evolution proceeds, thereby indicating the sense of time. However, 
the fact that such a relationship is no longer valid in the presence of local 
anisotropy of the pressure and/or dissipative processes and/or electric charge. 
This has already been discussed [321 ES] explaining its failure in scenarios 
where the above-mentioned factors are present [40] . Here we see how the 
electric charge distribution affects the link between the Weyl tensor and 
density inhomogeneity, suggesting that electric charge (whenever present) 
should enter into any definition of a gravitational arrow of time. 



3 Dynamics of Viscous Adiabatic Case 

In this case, heat flux vanishes, also, we assume that radiation density is 
zero and hence dissipation is only due to shearing viscosity. The energy- 
momentum tensor is obtained by replacing q a = e = in Eq. fl2.2p . Similarly, 
the Einstein-Maxwell field equations are found by using q = e = in the 
corresponding equations derived for non-adiabatic case. For junction condi- 
tions, the line element for the exterior region is taken as plane symmetric 
Reissner-Nordstrom spacetime given by 

— 2 /If p 2 —?M p 2 

ds\ = -{-f- + ^)dT 2 + (-^ + e -)~HZ- + 2?{dX* + dY% (3.77) 

where (x +0 , X +2 5 X +i ) — (T, X,Y, Z). The equation of hyper surf ace in 
terms of exterior coordinates is 

k + (T,Z) = Z-Zx(T) = 0. (3.78) 

Using Eq. fl3.78l) . the exterior metric over the hypersurface becomes 

,,,, r ,-2M p 2 . ,-2M e 2 . ,iZ SaiIm2 . „, 

+ Zl(dX 2 + dY 2 ). 
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The first junction condition yields the following equations 

dt 1 , 

Tr = A' Z - = B - < 3 ' 80) 



dT -2M e\i,,-2M e\, ,dZ^ s2 ,-, , , 

Tr = <^ + zI ), lhsr + zI ) "<lff ) |T - (3 - 81) 



Zjj, Zj£ ^ S 

Equation (I3.8ip implies that 



1 _2M e 2 

th* = NcH* - ±d%, N={^L + y. (3.82) 



The unit normal in terms of exterior coordinates is given by 

n+ = (-i E ,0,0,f), (3.83) 
The non-null components of the extrinsic curvature for the exterior spacetime 



are 

T ,+ _ l dZd 2 T d 2 ZdT N dN dT 3 3 dN dT dZ 2 

K+ = [ZN^-h = K+ 2 - (3.85) 
The second junction condition yields 

M = m(t,z) s = e. (3.86) 

The rate of change of mass with respect to T and Z are given by the following 
equations 

D f m(t, z) = -4ir(P z - ^ V F)UZ 2 + 5^ - 1), (3.87) 

2 

D~ z m(t, z) = An^Z 2 + Id 2 s + ^(/i 2 , - 1). (3.88) 

The description of these equations is the same as for the non-adiabatic case. 
Similarly, the dynamical equations can be obtained using q = e = in 
Eqs. fl232> f T23Tj) . For this case, Eq. (12371) becomes 

4 4 77? 4 ~ s 2 

(// + P Z --77F)D # C7 = -( /J + P 2 --r / F)[^ + 4vT(P 2 --7 / F)Z-^] 

-£ 2 [/^(P 2 -^F) + |(P z -P ± -2r 7 F)-^Sp]. (3.89) 
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The inequality given in Eq. ( 12.58]) remains the same. The interpretation of the 
terms in Eq.f l3.89]) is similar to that of the non-adiabatic case just excluding 
the effects of heat flux and radiation density. As there is no heat flux, so no 
transport equation is needed. The relationship between the Weyl tensor and 
energy homogeneity also remains the same. 

4 Conclusions 

Gravitational collapse is an outstanding phenomenon in gravitation theory. 
The aim of this work is to analyze the dynamics of gravitational collapse for 
plane symmetric configuration of real fluid. The conclusions are given in the 
following 

1. The junction conditions for both cases yield that masses of the interior 
and exterior regions are equal if and only if their corresponding charges 
are equal. For the non-adiabatic case, junction conditions also give 
a relationship between heat flux, effective pressure in z-direction and 
charge over the hyper surf ace ( I2.4ip . This equation implies that if s = 0, 
then the effective pressure in z-direction and heat flux are equal over 
the hypersurface. 

2. For the non-dissipative case, Eq.f l2^6]) reduces to £>x(§) = which 
shows that collapse is homologous, i.e., all the matter falls inward in a 
similar pattern. 

3. Condition for hydrostatic equilibrium is given by Eq.f l2.59j) . If this 
condition holds, then the collapsing process will stop and matter attains 
an equilibrium state. 

4. It is observed that charge will increase the rate of collapse if -| > D^s. 
Thus the chance of becoming a black plane increases in this case. 

5. In the non-adiabatic case, the radiation density increases inertial and 
active gravitational masses. Also, the outflow of radiation causes an 
increase in the rate of collapse and hence the collapsing process is ex- 
pected to be faster than the adiabatic case. 

6. The substitution of transport equation f|2.6ip in dynamical equation 
(I2.57P yields an additional factor a. This a term affects the inertial 
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mass density and gravitational force term. As a increases, the affected 
terms are decreased by the same amount and vice versa. The corre- 
sponding terms for a come from the term a^T in Eq.f l2.60j) which is 
the Tolman's inertial term. Hence the inertia of heat by increasing a, 
causes a decrease in inertial mass and gravitational force |JT]. Thus we 
can conclude that 

• If a — y 0, then inertial density and gravitational force are not 
affected by coupling. 

• If < a < 1, then inertia of heat causes a decrease in inertial and 
gravitational mass densities. 

• If a — > 1, then mass densities approach to zero. 

• If a > 1, then the gravitational force term becomes negative im- 
plying that reversal of collapse occurs due to the inertia of heat. 

7. Under certain conditions homogeneity in energy density and conformal 
flatness of spacetime are necessary and sufficient condition for each 
other. 

8. A relation (12. 75j) has been obtained exhibiting the way in which electric 
charge affects the link between the Weyl tensor and density inhomo- 
geneity. 
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